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Abstract 


Injuries  caused  to  the  head  by  ballistic  shock  waves 


during  bl 


acts  are  not  well  un¬ 


derstood.  It  has  been  postulated  that  traumatic  brain  injugy  (TBI)  can  be  caused  when  the 
blast  wave  causes  the  kinetic  energy  and  the  pressure  in  the  main  blood  vessels  to  oscillate 
rapidly  and  travel  to  the  brain  thereby  damagin g  The  Hkmal  fibers  and  neurons  or  with  the 
direct  wave  transmission  to  the  head.  In  the  direct  wave  transmission,  the  severity  of  blast 
wave  impact  can  be  reduced  by  ballistic  bAM^The  ballistic  helmets  currently  used  in  the 
military  have  been  designed  to  providmnroteCTion  against  penetrating  ballistic  projectiles  and 
their  effectiveness  against  blast  s]»kjrave  has  not  been  thoroughly  understood.  This  re¬ 
search  would  focus  on  developing  a  multiscale  computational  model  of  blast  impact  response 
of  high-performance  nap|^nposite  materials  for  the  helmet,  followed  by  estimation  of  blast 
energy  transfer  to  tMftissues  in  the  human  head.  Such  a  combined  atomistic-computational 
model  of  ballistic  response  of  nanocomposites,  coupled  with  a  human  morphology-specific 
computational  riw^rel  is  necessary  to  study  the  mechanics  of  blast  impact  on  human  head. 
Keywords:  kl  Hast  -induced  Traumatic  Brain  Injury,  Multiscale  Analysis,  Mathematical 
Homage  ruination,  Molecular  Dynamics 


The  skull  protects  the  brain,  which  is  the  most  sensitive  organ,  from  injuries  during 
fall,  stroke  or  blast  wave  impacts.  Traumatic  brain  injury  (TBI),  or  intra-cranial  injury,  is 
the  damage  caused  to  the  brain  by  external  mechanical  forces,  resulting  in  permanent  or 
temporary  impairment  of  the  brain  functions  [1].  Blast  injury  is  caused  by  a  sudden  increase 
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in  air  pressure  and  causes  injuries  in  spaces  containing  gas  [2],  In  a  blast,  shock  wave  travels 
at  the  speed  of  sound,  and  air  is  accelerated  by  this  shock  wave,  forming  a  high-velocity 
blast  wind.  It  has  been  found  that  during  Operation  Enduring  Freedom  and  Operation  Iraqi 
Freedom,  improvised  explosive  devices  (IEDs),  and  roadside  car  bombs  caused  aboq^teo 
of  US  casualties  in  Iraq  and  about  50%  in  Afghanistan  [3].  Blast- induced  now 

become  extremely  significant  as  studies  have  shown  that  between  January  2003  2005,  out  of 
the  450  soldiers  admitted  to  Walter  Reed  Army  Medical  Center,  59%  M^Eiagnosed  with 
TBI  of  which  56%  were  considered  moderate  to  severe,  and  44%  werd  considered  mild  [4,  5]. 

There  are  several  causes  of  blast-induced  TBI.  Primary  bias  t  injury  is  caused  by  a  high- 
pressure  shock  wave  interacting  with  the  body  [3] .  This  shigk  wave  produces  high  pressure, 
which  can  be  amplified  due  to  the  impedance  mismajmkJyween  the  skull  and  air  [6].  The 
direct  impact  of  the  shock  wave  causes  accelerate,  or  the  head,  and  since  the  inertia  of 
the  brain  is  less  than  the  inertia  of  the  skull,  the  brain  continues  to  vibrate  after  the  skull 
has  stopped  moving.  Cavitation  bubblei^prthen  formed  from  the  intracranial  pressure 
changes  [3]  and  these  bubbles  caqfe  il®ury  to  the  brain  as  they  collapse.  In  cases  when 
the  head  is  not  directly  subjogj^  tp*the  blast,  the  shock  wave  can  be  transferred  through 
the  abdomen  to  the  blopd  veSeis  or  other  fluid  pathways  in  the  thoracic  region.  This 
compression  effect  c^'Muse  waves,  which  transmit  kinetic  energy  to  the  brain  [7].  Another 
cause  of  traumj^^bram  injury  is  due  to  shrapnel  and  debris  from  explosions.  This  is  called 
secondary  MSsUihjtiry,  and  can  often  go  further  than  the  primary  high  pressure  blast  wave. 
The  arfe»n  of  the  different  body  parts  by  the  blast  wind  can  cause  tertiary  blast  injuries 
as  wUte®ue  to  the  difference  in  inertia  between  the  skull  and  brain,  higher  shearing  strains 
form  in  the  intra-cranial  region,  causing  shearing  injuries  [3].  Further,  the  high  pressure 
gases  following  the  explosion  are  at  very  high  temperatures  and  these  hot  gasses  can  cause 
quaternary  injuries,  which  include  pulmonary  injuries  from  the  toxic  gas  and  burns  and  the 
brain  may  also  be  damaged  if  the  skull  is  heated  excessively. 

Various  models  have  been  created  to  simulate  blast  waves  impacting  a  human  head.  A 
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finite  element  (FE)  study  was  initially  carried  out  that  could  be  parameterized  to  allow  for 
various  head  sizes  [8].  Grujicic  et  al.  demonstrated  the  use  of  a  Lagrangian-domain  model 
for  a  human  head  to  simulate  blast-induced  traumatic  brain  injury.  They  obtained  results 
for  pressure  and  stress  measured  at  different  points  in  the  head  model  [9].  Experimental 
tests  often  demonstrated  the  need  of  helmets  which  can  absorb  blast  pressurq^JT  thus 
mitigate  the  effects  of  blast  impact  on  the  human  brain.  At  the  same  tirn^ik^^nelmets 
must  protect  the  head  from  blunt  trauma  caused  by  shrapnel-induced  j^pUnds.  Ganpule  et 
al.  performed  finite  element  analysis  on  the  role  of  the  helmet  in  mitigating  blast  shock 
wave  propagation  [10].  They  concluded  that  the  highest  reflecde^i^pfrressure  when  a  shock 
wave  impacts  a  human  head  are  in  the  regions  of  concavity,  especially  at  the  nasion,  which 
is  the  nose-eye  cavity.  Kulkarni  et  al.  carried  out  a  cppipafative  study  on  ballistic  helmets 
[3].  This  review  included  helmet  materials  such  is  Kevlar  K29,  K129  fibers,  and  thermoset 
resins,  and  possible  future  materials  to  be  u^m^£o*helmets,  such  as  thermoplastic  polymers 
and  nano-composites. 

The  use  of  Kevlar-Carbon  Nap^^^based  composites  for  use  in  helmets  is  hypothetical; 
however,  there  have  been  nurner ohs  Jsfudies  on  Kevlar-nanocomposites  [11]  for  various  appli¬ 
cations  like  velocity  impact  studies  and  it  was  found  that  the  energy  absorption  and  impact 
response  characteris|rcS)f  the  composite  is  increased  due  to  the  presence  of  nanotubes  [12] . 
Similar  trends  «Br,g;alsw  observed  by  Randjbaran  et  al  [13]  for  the  ballistic  energy  absorption 
of  hybrid  composite'  laminates  reinforced  with  Kevlar  and  nanotubes. 

N a^focomposite  materials  are  often  better  alternatives  to  traditional  composite  materials 
becaulasm  their  excellent  mechanical  properties:  high  stiffness-to- weight,  strength-to-weight 
™pos,  and  excellent  structural  deformation  characteristics.  However,  the  major  deterrent 
to  experimental  investigation  of  novel  nanocomposites  in  the  helmet  design  is  the  extreme 
cost  and  the  lack  of  understanding  of  material  response  under  extreme  conditions.  This 
handicap  can  only  be  overcome  by  developing  multiscale  computational  models  that  consider 
the  influence  of  materials  at  the  nanoscales  and  macroscales.  In  this  work,  we  develop  a 
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multiscale  computational  model  of  the  helmet  material  with  reinforced  carbon  nanotubes. 
The  nanocomposite-based  helmet  material  would  be  studied  to  understand  the  effect  of 
energy  transfer  to  the  brain  under  blast  loading  conditions.  The  rest  of  the  paper  is  as 
follows,  Section  2  consists  of  the  development  of  the  computational  model  of  the  tywan 
head  and  the  representation  of  the  blast  wave.  This  section  is  followed  by  the  discussion  of 
numerical  results  on  the  amount  of  energy  transfer  from  the  blast  to  the  braip  in  a  protected 
and  unprotected  head  and  also  with  a  nanocomposite  helmet.  The  paj^r  concludes  with  a 
discussion  of  results. 


pamr  c 


2.  Ideal-blast  Waves 


There  are  many  sources  of  blasts  that  occur  in  a  battle-field  that  produces  pressure  waves 
by  rapid  release  of  energy.  Explosion  sources  alsowMisduring  the  rapid  vaporization  of  thin 
metal  films  or  fine  wires.  However,  more  well-known  blast  sources  are  from  chemical  or  nu¬ 
clear  reactions,  when  a  sudden  increase  in^%ir  pressure  causes  a  shock  wave  to  form.  This 
shock  front  is  almost  instantaneous  and  causes  nearly  discontinuous  increases  in  tempera¬ 
ture,  density,  and  pressure  [14]'  A  ^pressure  wave  can  be  approximated  using  a  Friedlander 
equation,  which  can  be  y^»<to  describe  the  physical  properties  of  an  ideal  blast  shock  wave 
as  given  by  Equation  (1)  [15,  16]  (see  Figure  1).  The  most  efficient  means  of  representing 
a  blast  impact  is  by  The  use  of  data  from  ConWep.  ConWep  is  a  collection  of  calculation 
routines  tl^Tcharacterizes  the  effect  of  conventional  weapons  that  includes  air  blast  routines, 
breaqj^^rj^rmg,  ground  shock,  and  fragment  and  projectile  penetrations  etc.  [17,  18]. 


rP 


(i) 


3.  Computational  Model  of  Human  Head  and  Helmet 

The  3D  FE  computational  models  are  developed  from  the  medical  image  data  to  create 
STereoLithography  (.stl)  files  [19].  The  FE  model  of  the  human  head  was  obtained  from  the 


4 


ACCEPTED  MANUSCRIPT 


MR/CT  data  obtained  from  the  Visible  Human  Female  database  from  the  National  Library 
of  Medicine  [20].  This  study  partially  follows  the  process  outlined  by  Vonach  et  al.  for 
creating  the  model  of  the  human  head  [19].  Following  the  creation  of  the  geometric  model, 
tetrahedral  or  hexahedral  meshing  algorithms  are  used  to  generate  the  equivalent^^^te 
element  mesh.  The  MR/CT  images  were  then  converted  to  a  FE  mesh  using  segmentation 
followed  by  smoothening  of  the  resulting  STL  files  and  final  conversion  into  ,  Ih.  FE 

analysis  was  carried  out  using  the  general  purpose  FE  code  Abaqus  (/j^aMsVersion  6.12, 
Dassault  Systmes)  to  study  the  effects  of  blast  induced  pressure  differences.^  the  human  head 
with  and  without  an  helmet,  using  ConWep  blast  loading.  TheCk>h]pep  implementation  in 
Abaqus  is  used  for  the  simulations.  The  FE  geometry  ofi|he  human  head  model  includes 
the  brain,  surrounding  cerebral  spinal  fluid,  and  the  skull.  In  the  FE  model,  the  head  is 
fixed  at  the  base  of  the  neck  with  encastre  boundary  conditions.  The  model  consists  of 
the  skull  having  5998  elements,  the  cerebral  ipinal  fluid  with  5352  elements,  and  the  brain 
model  contains  6745  elements,  for  a  total  ot|^o95  elements  of  type  C3D4.  In  this  study,  the 
heterogeneity  of  the  brain  materiaj^rmms  non-isotropic  behavior  is  not  taken  into  account. 
A  Mooney-Rivlin  model  is  for  the  brain  [21]  and  the  strain  energy  function  is 

given  by  [22,  23]: 


—  Cio (Ji  ~  3)  +  C0i(J2  —  3)  +  — ( Jei  —  l)2 


(2) 


where  W  is  the  Strain  energy  potential,  Jej  is  the  elastic  volume  ratio,  J\  is  the  first  invariant 
of  the  de-vftitoric  strain,  J2  is  the  second  invariant  of  the  deviatoric  strain,  and  Ci0,  Coi,  and 
material  constants  dependent  on  temperature.  Linear  viscoelasticity  [24]  is  used  to 
estimate  the  second  Piola-Kirchhoff  stress  by  the  following  convolution  integral  [21]: 


Sij  J  Cljkl{t-r)5^ldr  (3) 

o 

where  is  the  second  Piola-Kirchhoff  stress,  Gijki(t  —  r)  is  the  relaxation  modulus  for  the 
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different  stress  measures,  and  Ekl  is  the  Greens  strain  tensor.  The  Cauchy  stress  cr  (c %)  can 
be  obtained  from  the  Second  Piola-Kirchoff  stress  from  [22] 


r  =  -F  •  S  •  Ft 


where  F  is  the  deformation  gradient  tensor,  and  J  is  the  Jacobian  of  transformalji^^^^nally, 
this  stress  is  superimposed  to  the  stress  tensor  which  was  determined  from  the_stram  energy 
function.  A  Prony  series  can  be  used  to  represent  the  relaxation  modulus  ror  an  isotropic 
material  as  follows  [21]: 

G{t)  =  G0  +  J2  Gie~Pit  (5) 


where  /3j  is  the  decay  constant  and  Gt  is  the  relaxation  modulus.  The  parameters  selected 
for  the  brain  are  density  of  1040  kg/m 3  [9]  witte^ooney-Rivlin  coefficients  C\0  =  3102.5 
Pa,  C0 1  =  3447.2  Pa,  and  D%  =  1.095  CPa^r* bulk  modulus  of  2.19  GPa  [21],  The 
cerebral  spinal  fluid  is  modeled  using  flurefiple  properties  using  the  linear  Us-Up  equation 
of  state  (EOS),  with  the  final  Hugoniot(  pressure  used  as  the  constitutive  relationship.  Us 
is  the  shock  velocity  and  Up  jdrzke  jmrticle  velocity  and  it  can  be  considered  that  a  linear 
relationship  exists  between  the  Mock  velocity  and  particle  velocity  given  by  Us  =  c0  +  sUp. 
In  this  linear  relatioj^ftr,  p  =  — Ki^ ,  where  K  is  the  bulk  modulus  of  elasticity,  p  is  the 
pressure  rate,  an m  iij  #  the  deviatoric  strain  rate  [21],  The  linear  form  is  given  by: 


(1  -  srj)  2 


(6) 


linear  response  described  here  is  obtained  by  setting  the  parameters  s  =  0.0  and 
To  =  0.0  since  K  =  p0cl  for  small  strains.  This  gives  a  c0  of  1451.127  m/s  using  a  bulk 
modulus  of  2.19  GPa  [21]  and  the  density  is  taken  as  1000  kg/m3.  The  skull  is  considered  as  a 
linear  elastic  material  and  is  modeled  based  on  Mori-Tanaka(M-T)  [25]  based  homogenization 
of  a  porous  bone  structure.  For  an  elastic  composite  material,  the  effective  constitutive 
relations  are  given  by  the  volume  average  of  the  stress  and  strain.  Similarly,  for  each  phase 
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on  the  microscale,  the  constitutive  relation  can  be  given  as: 


(°)t  =  Ck(e)t 

where  {.)^ot  is  the  volume  averaged  state  of  phase,  including  the  matrix,  fiber 
interphase  layers  [26].  The  volume  averaging  (.)  of  the  state  variables  for 
continuum  (EC)  are  given  by 


(°*t 


--  ~  J  vapdv-  {eap)EC=  y  J  eaPdv; 


(7) 

id  any 


.ivalent 


(8) 


and  for  an  N- particle  atomic  ensemble  we  have:  (eap)  —  j?  Z)  £ad- 

The  porosity  of  the  skull  bone  is  taken  as  0.8Tv®fh  is  that  of  a  cancellous  bone  and 
having  a  Poisson’s  ratio  of  0.32,  with  a  Young’s  modulus  of  22.0  GPa  [27]  for  the  matrix 
material.  The  density  of  the  skull  bone  iijrai^i  as  1412  kg/m3  [9].  The  effective  Young’s 
modulus  and  Poisson’s  ratio  of  the  co^fcsite  bone  (porous)  tissue  using  Mori- Tanaka  method 
is  given  by  Equations  (9)-(12^f^8]^^a  the  homogenized  modulus  of  the  skull  is  taken  as 
2.174  GPa,  with  Poisson’§jfetfjf|glT.253  and  porosity  of  0.82  [27]. 


where 

where,  7o  -  ^±4. 


.4 


E  —  2p\l  + 


3K  -  2JI  1 
2(3K  +  jl)\ 


:  =  if0{i 
n  =  no  |i 


cjKj-Ko) 
"  K0  +  37o(l  —  c)(K 

,  c(ni  -  no) 

no  +  2^o(i -  c)(n  1 

_ B 

~  3(1- 


'  Ahi)  J 


(9) 


(10) 


(11) 


3  K-E 


(12) 
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The  helmet  made  of  thermoset  resin/KevlarK129  composite,  with  density  1440  kg/m3, 
effective  Young’s  modulus  of  1.24  GPa,  and  Poisson’s  ratio  of  0.36  is  used  in  the  simula¬ 
tions  [29]  (as  shown  in  Table  1).  The  inner  padding  is  considered  to  have  a  density  of 
136%/m3,  with  a  Young’s  modulus  of  8.0  MPa,  and  Poisson’s  ratio  of  0.2  [29],  Additiphstlly 
we  consider  a  carbon  nanotube  reinforced  thermoset  resin  also  to  understand  j^ejftiergy 
dissipation  characteristics  of  the  helmet  material. 

Table  1:  Properties  of  skull,  brain  and  helmet  material  used  in  computa^piSfcgimulations 


Material 

Elastic  (E)/Bulk  (K) 
Modulus 

Poissonls'd 

ratio 

iy^isity 

1  (kg/m3) 

Brain  Tissue 

K  =  2.19  GPa 

W 

1040 

Skull 

E  =  2.174  GPa 

0.2b3 

1412 

Helmet  (KevlarK129) 

E  =  1.24  GPa 

iC36 

1440 

The  Young’s  modulus  of  a  (5,  5)  single- wallmj^bon  nanotube  was  calculated  using 
LAMMPS,  a  molecular  dynamics  simulator-  It  should  be  emphasized  that  atomistic  sim¬ 
ulations  are  carried  out  only  on  the  CNT%nd  not  on  a  Kevlar-CNT  ensemble  and  it  is 
assumed  that  the  properties  of  the  CNT  obtained  from  the  MD  simulations  would  bridge 
the  scales  from  the  atomistic  to  the  continuum.  In  MD  simulations,  the  motion  of  particles 
in  an  atomistic  system  irg«erned  by  the  Hamiltonian,  which  is  a  function  of  the  position 
and  momentum  of  ttejrarticles,  and  the  Hamiltonian  equations  of  motion.  The  Hamiltonian 
represents  the  mtal  energy  of  an  isolated  system  and  is  composed  of  the  sum  of  the  potential 
and  kinetic  enefigy  and  the  thermodynamic  terms 


(13) 


where  U  (rN)  is  the  potential  energy  from  intermolecular  interactions  as  a  function  of  the 
spatial  ordinate  rN ,  A-  pf  represents  the  momentum  p  of  the  particle  i  with  mass  rnt ,  which 
is  a  function  of  the  absolute  temperature  called  the  kinetic  energy.  The  time  derivative  of 
the  Hamiltonian  gives 
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dR  1  dU 

—rr  =  —  >  Pi-Pi  +  >  ^—-n  =  0 

dt  m  dri 

and  with  the  spatial  derivative  we  obtain  the  equation  of  motion 


dR 

d/rt 


dU 

dri 


(14) 


(15) 


)f  valem 


.bond  energies 

(16) 


The  total  potential  of  the  nanotube  structure  is  given  by  the  sum  of 
and  the  nonbonding  interactions  [30,  31,  32] 

»“=EEh”+^V 
3  3>l 

where  V®  is  the  potential  energy  due  to  bondinsym^^-^  is  the  potential  energy  due  to 
nonbonding  interactions.  The  reliability  of  an  MD  si  mulation  depends  mainly  on  the  type  of 
potential  functions.  Empirical  bond-orde^pi^nuals  of  the  Reactive  Empirical  Bond  Order 
(REBO)  Potential  of  Brenner  type^S^found  to  be  excellent  for  the  simulation  of  carbon 
based  materials  [33,  34,  35] .  Th^se  l§|ires  of  potentials  are  capable  of  describing  the  changes 
in  the  bonding  in  a  systena^but  lacks  in  describing  long  range  interactions,  which  are  often 
neglected  [36].  The  Adaptive  Intermolecular  Reactive  Empirical  Bond  Order  (AIREBO) 
Potential  has  b<pi  used  in  the  simulations  [37] .  The  total  potential  energy  of  an  atom  in 
an  AIREBQapotentfal  is  given  by 


IBQjaofenth 

'•  to 

S^E.' 


i+ev+y:  e  e™?s,on 


(17) 


fiREBO _ 

k/hj 

I  EBEBO  represents  the  Reactive  Empirical  Bond-Order  (REBO)  Potential  and  E[d  rep¬ 
resents  the  non-bonded  interactions  and  E^lRSION  represents  the  torsional  component  of  the 
potential  function.  The  force  of  attraction  and  repulsion  (Faf})  experienced  by  each  molecule 
is  obtained  from  the  gradient  of  the  potential  field. 
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f„r  =  - 


dUtot 

dra8 


(18) 


The  molecular  dynamics  simulations  included  a  minimization  of  the  input  systerri^hn 
equilibration  to  300  K  for  30000  steps,  and  a  deformation  test  for  20000  steps  at  a  stra% 
rate  of  0.01  ps_1.  Each  time  step  was  taken  as  0.2  fs  and  uses  a  reactive  A I REBO -potential 
function  developed  by  Stuart  et  al.  for  the  modeling  of  intermoleculai^Mg^fions  and 
chemical  reactions  in  a  condensed-phase  hydrocarbon  system  [37].  Yrong’s  modulus 
of  the  carbon  nantoube  was  estimated  to  be  1.004  TPa  and  it  .closely  matches  the  values 

J3 

reported  in  literature  of  1  TPa,  with  a  Poisson’s  ratio  of  0.14  [38]  and  having  a  density  of 
1330  kg/m 3  [39].  A  5%  volume  fraction  of  the  carbon  nanotube  was  assumed  in  the  analysis 
and  the  resulting  composite  is  considered  to  have  an  effective  Young’s  modulus  of  1.3772 
GPa,  an  effective  Poisson’s  ratio  of  0.35,  withm  cnrisity  of  1434.5  kg/m3. 


4.  Numerical  Results  and  Discussion 


Three  different  magnitudes^:  explosion  are  modeled  using  the  ConWep  blast  loading 
pressure  wave  as  shown  in  Figure  z,  at  a  distance  of  0.8  m  in  front  of  the  head  model.  The 
coup  side  correspond^ to  tnwsite  of  impact  of  the  explosion,  which  in  this  case  is  the  front 
of  the  brain.  + 

The  maxmilim  shear  stress  was  found  to  be  23.2  kPa  for  0.038  kg  blast  loading  and 
with  0.093  kg  of  TNT,  the  maximum  shear  stress  was  obtained  as  53.4  kPa.  For  0.227  kg 
explosion,  the  maximum  shear  stress  was  obtained  as  90.9  kPa  and  these  compares  well  with 


wa&pinulations  carried  out  by  Chafi  et  al.  [21],  with  increasing  discrepancy  for  larger  blast 
kmding,  and  are  attributed  to  the  difference  in  the  model  development  from  the  MRI  data. 
The  coup  shear  stress  for  the  brain  for  the  three  different  explosions,  measured  with  respect 
to  time  is  indicated  in  Figure  3.  Simulations  are  further  carried  out  for  the  same  magnitudes 
of  TNT  exploding,  at  a  distance  of  0.8  m  to  the  rear  of  the  head.  The  roles  of  the  coup  and 
countercoup  sites  are  now  reversed  as  the  detonation  is  now  performed  at  the  rear  of  the 
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head  and  as  shown  in  Figures  4,  the  maximum  principal  stress  at  the  coup  and  countercoup 
sites  are  found  to  be  significantly  higher  at  the  coup  site. 

A  helmet  is  now  added  to  the  head  model  with  a  padding  on  the  inside.  Similar  simula¬ 
tions  were  carried  out  to  demonstrate  the  efficacy  of  the  helmet  in  mitigating  the  effedf%.)f 
shock  waves  impacting  the  human  head.  The  results  obtained  from  detonations  omwng  at 
the  front  of  the  head  are  compared  to  an  unprotected  head  model  and  the  diffe'rei^trnagni- 
tudes  of  the  maximum  principal  stress  were  obtained  for  pristine  helmuLinaterial  as  shown 
in  Figure  5.  The  principal  stresses  are  compared  from  an  unprotected  hea^r  to  that  of  a  pro¬ 
tected  head  as  shown  in  Figure  6  and  it  can  be  seen  that  for  highest  load  of  0.227  kg  TNT, 
the  maximum  principal  stress  is  significantly  reduced  for  thp, protected  head  when  compared 
to  the  unprotected  head.  This  difference  was  not  significant  for  higher  loads,  which  indicate 
that  greater  amount  of  stresses  are  transferred  tf&tehe  brain  tissue  as  the  impact  pressure 
gets  higher. 


The  efficacy  of  the  helmet  in  protecting^^brain  from  blast- induced  TBI  can  be  measured 
in  terms  of  the  transfer  of  strain  ra^rg^density  (SED).  SED  is  the  amount  of  strain  energy 
stored  per  unit  volume.  Simulations  were  performed  to  calculate  the  energy  density  for  the 
three  different  magnitudes  of  explosions.  Figures  7,  8  and  9  indicate  the  energy  density  at 
the  front  of  the  brain  f^up  site)  for  the  porous  unprotected  head  model  (M-T  model  skull), 
the  head  protected  by  tin1  ACH,  and  the  head  protected  by  the  CNT-enhanced  helmet. 
It  can  be  Jurat  the  energy  density  was  greatly  reduced  when  the  unprotected  head  is 
compared  to  the  head  with  the  advanced  combat  helmet  and  the  carbon  nanotube-enhanced 
helmet  material.  This  indicates  that  significant  energy  is  absorbed  by  the  helmet  and  less 
s|pck  energy  per  unit  volume  reaches  the  brain.  However,  the  same  discrepancies  discussed 
previously  are  seen  when  comparing  the  CNT-enhanced  helmet  with  the  original  Advanced 
Combat  Helmet  (ACH).  This  discrepancy  can  be  rectified  by  considering  the  use  of  an  energy 
dissipation-based  material  model  for  the  analysis. 
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5.  Conclusions 

Head  injuries  by  ballistic  shock  waves  during  blast  impacts  are  not  well  understood  and 
the  damage  to  the  brain  tissue  leads  '  ’ 


slow  to  appear.  The  ballistic  helmets  currently  used  in  the  military  have  been  desk 


protection  against  penetrating  ballistic  projectiles  and  their  effectiveness  again^ 
has  not  been  thoroughly  understood.  The  major  outcome  of  this  study  isJii^ck 
of  an  enhanced  computational  model  that  helps  in  our  understanding  of  blast-in 


and  exploring  the  possibihty  of  mitigation  of  head  injury  using  mjenocomposite 


helmets. 
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Figure  1:  Representation  of  Friedlander  Blast  Wave 


,(j 


Figure  2:  Conwep  pressure  applied  at  front  of  skull 
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Figure  4r  Principal  stress  in  the  brain  at  the  (a)  coup  site  and  (b)  counter  coup  site  -  for  an 
m^^ected  rear  detonations 
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Figure  6:  Maxppra®  principal  stress  with  and  without  helmet  protection  in  the  brain 
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